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ABSTRACT 

We present HST/NICMOS photometry, and low-resolution K-band spectra of the GLIMPSE9 stellar 
cluster. The newly obtained color-magnitude diagram shows a cluster sequence with H-K^ =^ 1 
mag, indicating an interstellar extinction ^Ka= 1.6±0.2 mag. The spectra of the three brightest stars 
show deep CO band-heads, which indicate red supergiants with spectral type M1-M2. Two 09-B2 
supergiants are also identified, which yield a spectrophotometric distance of 4.2 ± 0.4 kpc. Presuming 
that the population is coeval, we derive an age between 15 and 27 Myr, and a total cluster mass 
of 1600 ± 400 Mq, integrated down to 1 Mq. In the vicinity of GLIMPSE9 are several HIT regions 
and SNRs, all of which (including GLIMPSE 9) are probably associated with a giant molecular cloud 
(CMC) in the inner galaxy. GLIMPSE9 probably represents one episode of massive star formation in 
this CMC. We have identified several other candidate stellar clusters of the same complex. 

Subject headings: stars: evolution — infrared: stars 



1. INTRODUCTION 

An understanding of the mechanisms of formation and 
evolution of massive stars is of broad astronomical inter- 
est. Through mass- loss and supernova explosions, mas- 
sive stars return a significant fraction of their masses to 
the interstellar medium (ISM), thereby chemically en- 
riching and shaping the ISM. Being very luminous, they 
can be identified in external galaxies, providing spec- 
trophotometric distance. Observational constraints on 
the formation and evolution of massive stars are, how- 
ever, difhcult to obtain due to the rarity of these objects 
and their location in the Galactic plane, where interstel- 
lar extinction can hamper their detection. 

Massive stars can be identified by their ionizing radi- 
ation, which creates easily identifiable HII regions. Fur- 
thermor e, since the ma jority of massive stars are born in 
clusters (jLada fc Ladai r2003). they are also identified by 
locating young massive stellar clusters. Over the past 
decade, infrared and radio observations of the Galac- 
tic plane have revealed several hundred new HII regions, 
more than 50 new candidate supernova r emnants (SNR) 
. (iGiveon et al.l 120051: iHelfand et al.ll200l. and 1500 new 
candidate infrared stellar clusters fe.g. iBica et al.l [20031 : 
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iMercer et all 120051 : iFroebrich et al.ll2007t ). which are of- 
ten found in the direction of HII regions. Only a few 
of these infrared candidate clusters have been confirmed 
with spectro-photometric studies; the analysis is often 
restricted only to stellar clusters and does not include 
the cluster environment. A combined study of stellar 
clusters and their associated molecular clouds is a pow- 
erful tool to understand star for mation. Clusters a ppear 
to form in large complexes (e.g. ISmith et al.ll2009r ). The 
temporal and spati al distribution of clusters varies from 
cloud to cloud (e.g . iHomeier &: Alve£ll2005l : iKumar et al.l 
l2004t IClark et alTl2009a^ ■ indicating that external and 
internal triggers are both at work. Supernova explosions 
may trigger subsequent episodes of star formation in the 
same cloud. The presence of SNRs indicates that a cloud 
has already undergone massive star formation, and the 
study of stellar clusters associated with SNRs can shed 
light on the initial masses of the supernova progenitors, 
and therefore on the fate of massive stars. 

By locating new HII regions, and young stellar clus- 
ters, we also obtain information on large scale Galactic 
structures. So far, we have only located clusters that 
reside in the near s ide of the Galaxy, with a few ex- 
ceptions, e.g. W49 (|Homeier fc Alved [20051) . Many is- 
sues on Galactic structures are still open, e.g. the ex- 
act number of spiral arms, the lack of star formation in 
the central 3 kpc, and the possible existence of a ring 
of massive star formation surrounding the central bar. 
Stellar clusters selected from infrared observations are 
promising tracers for Galactic studies because they sam- 
ple a larger port ion of the Galactic plane than those from 
optical surveys ("Messineo et al."2009^ iDavies et alll2007l : 
.Figer et al. 2006: Clark et al. 2009b). 

The candidate cluster number 9 in the list by 
IMercer et all (|2005D (hereafter, GLIMPSE9) is an ideal 
target to study the issues above mentioned, because it is 
located in projection on a CMC that hosts several HII 
regions, supernova remnants, and candidate clusters. So 
far, we know only of one other Gala ctic cluster associated 
with a SNR ([Messineo et all 120081 ). This complex is on 
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Fig. 1. — The greyscale above shows a 2MASS /fg-band image, 
while contours show 24 fim emission from MIPSGAL. Contour lev- 
els are 100, 200, 300, 400, and 500 MJy/sr. The two boxes indicate 
the location and size of the two NIC3 fields. Extended emission at 
a 5 o" level is seen towards the field. 

the Galactic plane, at an heliocentric dista nce of ^ 4.9 
kpc ([Albert et all 120061 : iLeahv et al.]l2008l ). and longi- 
tude 1=~ 23°, and represents an episode of massive star 
formation in the direction of the inner Galaxy. Here we 
present a spectro-photometric study of the GLIMPSE9 
cluster. 

In Sect, m we present the available photometric and 
spectroscopic data, and the process of data reduc- 
tion. The spectral analysis and color-magnitude diagram 
study are given in Sect. [31 Information about the parent 
GMC and other candidate stellar clusters are presented 
in Sect. [H Finally, in Sect. [5] we summarize the results 
of our investigation. 

2. OBSERVATION AND DATA REDUCTION 
2.1. NICMOS/HST observations 

Images were taken with NICMOS on board the HST 
on July 9, 2008 as part of the GO program 11545 (P.I. 
Ben Davies). Two fields were observed; one centered 
on the GLIMPSE9 cluster (RA= 18l'34™09.89", DEC=- 
09?14f04.8") and another located at (RA=18l'34™12.76", 
DEC=-09?12.'45.4"). The latter, which we caU the con- 
trol field, was imaged in order to study the background 
and foreground stellar population in the direction of the 
cluster. The location of the two fields is shown in Fig. [TJ 
The NIC3 field of view (51.5"x 51.5") covers the cluster 
to a radius roughly equal to two half light radii (~30", 
as measured on 2MASS images). 

We used the NIC3 camera with the F160W and F222M 
broadband filters, and with the F187N and F190N nar- 
rowband filters. The fields were dithered by 5.07" in a 
spiral pattern (6 positions). The STEP2 sequence of the 
MULTIACCUM readout mode with 13 reads was used 
for exposures with the F160W filter, giving an integra- 
tion time of 19.94 s per exposure; the STEP8 sequence 
with 12 reads was used for exposures with the F222M fil- 
ter (55.94 s), and the STEP8 with 10 reads for exposures 
with the F187N and F190N filters (39.95 s). 
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Fig. 2.— Top: F160W mosaic of the GLIMPSE9 cluster. 
Squares and identification numbers (from Table [TJ show the lo- 
cation of spectroscopically observed stars. Stars #10 and #11 fall 
outside the field of view. Bottom: F160W mosaic of the control 
field. 

2.2. NICMOS/HST data reduction 

The images were bias subtracted, dark corrected, and 
fiat-fielded by the standard NICMOS pipefine CALNICA 
(see the NICMOS Data Handbook v7.0^). The six 
dithered exposures of each observation were re-sampled 
into a final mosaic with a pixel scale of 0.066". The clus- 
ter mosaics are shown in Fig. [2l 

A major photometric uncertainty is inherent for im- 
ages with the NIC3 camera, and is due to the combina- 
tion of an under-sampled point spread function (PSF) 
with the lower intrapixel sensitivity of the camera (up 
to 10-20%, see the NICMOS Data Handbook v7.0). The 
dithering observing strategy reduces this uncertainty. A 
photometric analysis of the mosaics was carried out with 
DAOPHOT Stetson (1987) within the Image Reduction 
and Analysis Facility (IRAF). We used the images of 
the control field, which is less crowded, to build a point- 
spread- function (PSF); seven isolated and bright stars 
were selected in the F160W and F222M images. Due 
to the small number of isolated stars no attempt was 
made to model a spatially varying PSF. Aperture pho- 
tometry with a radius of 1.1" was performed on bright 
and isolated stars, and the average difference between 
these magnitudes and the PSF-fitting magnitudes from 
DAOPHOT was measured. This aperture correction was 
applied to the whole catalogue. We obtained calibrated 
magnitudes in the Vega system with the transformation 

® |http://www.stsci.edu/hst/nicmos7] 



The GLIMPSE9 stellar cluster. 
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Fig. 3. — Magnitude errors from DAOPHOT for the cluster field 
([F160W] in the upper panel, and [F222M] in the lower panel). 

equations between counts/s to Vega magnitudes as given 
in the NICMOS Data Handbook. Astrometry and pho- 
tometric calibrations into the 2MASS system were ob- 
tained using a set of 13 point sources \ yith good qual- 
ity H and Ks measurements in 2MASS (jSkrutskie et al.1 
l2006f ). The reference stars span a, H — Ks range from 0.2 
to 1.4 mag, and a Ks range from 6.3 to 14.0 mag. 
The transformation equations are: 



Ks - 

(-0.01 ±0.10) + (-0.03 ±0.1 



[F222M] = 

8) X ([F16017] - [F222M]), 



H - [Fimw] = 

(-0.08 ± 0.09) ± (-0.22 ± 0.07) x ([F160W^] - [F222M]), 

with a standard deviation of 0.11, and 0.10 mag, re- 
spectively. [i^222M] and [F160M^] are magnitudes in 
the VEGA photometric system (following the NICMOS 
manual), and Ks and H are magnitudes in the 2 MASS 
system. A plot of the photometric errors given by 
DAOPHOT is shown in Fig. [31 

We also extracted point sources from the F187N and 
F190N images. The difference of [F187N] and [F190N] 
magnitudes show a scatter of 0.3 mag. Within this un- 
certainty, no emission lines were detected. 

2.3. Spectroscopic data 

We obtained spectroscopic observations with the Infra- 
Red Multi-Object Spectrograph (IRMOS) at the Kitt 
Peak Mayall 4m teles cope on September 25th, 2007 
(jMacKentv et all 120031) . We used the KIOOO grating in 
combination with the K filter to cover the wavelength 
region from 1.95 /im to 2.4 /^m with a resolution of 
R=~ 1000. A number of 28 exposures of 1 minute each 
were taken in two nodded positions. In order to remove 
variable background signals each exposure was followed 
by a dark observation of equal integration time. Neon 
lamp and continuum lamp observations were taken soon 
after the target observations. Dark-subtracted science 
frames were combined and flat-fielded. A two dimen- 
sional de- warping procedure was then used to straighten 
the stellar traces before extraction. Wavelength calibra- 
tion was obtained by u sing both neon lines and OH lines 
iva fc O riglia 1992^ A total of seven spectra were ex- 
tracted. Each target spectrum was divided by the spec- 
trum of an A2V star in order to correct for atmospheric 
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Fig. 4. — Left panel K-band spectra taken with NIRSPEC and 
IRMOS. The number labels refer to Table[T] The spectra have been 
corrected for interstellar extinction using Aj^^=1.5 mag and the 
extinction law by Mcssineo et al. (200^). For tree stars (#1, #4, 
and #5) with both KECK and IRMOS observations both spectra 
are shown (the Keck spectrum is the top one). 



absorption and instrumental response. The Br7 of the 
telluric spectrum was eliminated with linear interpola- 
tion. The intrinsic shape of the telluric spectrum was 
re moved by dividing it by a black body function of 9120 
K (|Bhim et al.ll2n0f)h . 

Additional spectroscopic observations were carried out 
with NIRSPEC at the Keckll telescope under program 
U050NS (P.I. M. Rich), on July 11th, 2008. We used the 



K filter and a 42" x 0.570" slit. A wavelength coverage 
from 2.02 /im to 2.45 /im and a resolution R=1700 were 
obtained. For each star, two exposures were taken of 10s 
each, in two nodded positions along the slit. We used a 
continuum lamp observation as a flat field, and Ar, Ne 
and Kr lamps observations for the wavelength calibra- 
tion. Pairs of nodded positions were subtracted and flat- 
fielded. Atmospheric absorption and instrumental re- 
sponse were removed by dividing each extracted target by 
the spectrum of a B0.5V telluric standard (HD1762488), 
and multiplying fo r a blackbody spectrum of 32,060 K 
(|Blum et al il2000f ). A total of seven spectra were ex- 
tracted from the NIRSPEC observations. Three of the 
stars observed with IRMOS were also observed with NIR- 
SPEC. 

A chart of the observed stars is shown in Fig. [21 The 
seven spectra with NIRSPEC and the additional four 
spectra with IRMOS are shown in Fig. [H In Table [H 
we list the coordinates of the stars spectroscopically ob- 
served. 



3. ANALYSIS 
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TABLE 1 

List of stars spectroscopically observed. 



ID 


RA 


DEC 


frame 


Spec, type 


EW(CO) 


H 






1 


18 34 09.266 


-09 14 00.74 


NIRSPEC(118r)/IRMOS(ir3) 


Ml-OI 


-43-39 


07.80 


07.17 


0.6 ±0.1 


2 


18 34 10.064 


-09 13 57.74 


NIRSPEC(1181) 


<G 




10.30 


09.20 


1.6 ±0.1 


3 


18 34 08.692 


-09 14 11.07 


NIRSPEC(1201) 


OBI 




08.93 


07.96 


1.3 ±0.1 


4 


18 34 08.537 


-09 14 11.83 


NIRSPEC(120r)/IRMOS(ir4) 


OB 




10.21 


09.14 


1.5 ±0.1 


5 


18 34 09.857 


-09 14 23.28 


NIRSPEC{122) /IRMOS(ir2) 


M1-2I 


-43-49 


08.43 


07.05 


1.7± 0.1 


6 


18 34 10.353 


-09 13 48.99 


NIRSPEC{124r) 


M5III 


-32 


10.41 


09.44 


1.0 ±0.1 


7 


18 34 11.348 


-09 13 46.47 


NIRSPEC{1241) 


<G 




12.05 


10.93 


1.6 ±0.1 


8 


18 34 10.352 


-09 13 52.95 


IRMOS(irl) 


M2.5I 


-49 


07.58 


06.30 


1.5±0.1 


9 


18 34 08.228 


-09 14 03.25 


IRM0S(ir5) 


<G 




10.77 


09.78 


1.4 ±0.1 




18 34 06.662 


-09 14 47.40 


IRM0S(ir6) 


>K0 




12.26 


10.26 


2.7± 0.1 


11=^ 


18 34 05.345 


-09 14 24.01 


IRM0S(ir7) 


<G 




10.05 


8.96 


1.6 ±0.1 



Note. — For each star, number designations and coordinates (J2000) are followed by the instrument name (plus frame 
name), spectral classification, EW(CO), H and Kg magnitudes obtained from the HST images, and the estimated interstellar 
extinction (in magnitude). 

2MASS magnitudes are listed for stars N. 10 and 11, because they arc outside the area covered by the NIC3 mosaic. 



3.1. Spectral types 

Spectral classification was perform ed by comparing; the 
spectra with spectral atlases (e.R. iHansoii et alj [l996l 
2005Hlva nov et al.ll2004l: I Alvarez et al.ll2000l: iBlum et all 
1996; Kl einmann fc Halllll98(il : [Wallace fc Hinkld 119961 ). 

X-band observations enabled us to classify both late- 
and early-type stars. Typically, spectra of late-type stars 
show CO band- head at 2.29 /im and atomic lines from 
Mg I, Ca I, and Na I; early-type stars can be identified 
by detecting hydrogen (H) lines, helium (He) lines, and 
other atomic lines, e.g. CIV triplet at 2.069, 2.078, and 
2.083 /im, and broad emission at 2.116 /im, which is due 
to NIII, CIII and Hel emission. 

Stars #1, #5, #6, #8, and #10 show CO bands in 
absorption, which indicate low effective temperatures. 
Since the absorption strength of CO band-heads in- 
creases with decreasing effective temperature Toff, but 
with increasing luminosity L, giant stars and supergiant 
stars follow different equiv alent width EW(C O) versus 
temperature relations (e.g. iDavies et al.ll2007l ). For each 
star, we measured the EW(CO) band-head feature be- 
tween 2.285 /im and 2.315 //m, with an adjacent con- 
tinuum measurement made at 2.28-2.29 /tm. Then, we 
compared these measurem ents to those of template stars 
einmann &: Halll 119861 ). and determined the spectral 
type (see Fig. [5]). Stars #1, #5 and #8 appear to be 
red supergiant stars (RSGs) with spectral types between 
MO and M2. Star #6 is most likely a giant M5 because 
of its small EW(CO) and of its lower luminosity {Ks = 
9.4 mag). Star #10 also shows CO band-head in ab- 
sorption, but because of poor signal-to-noise 10) no 
spectral typing is attempted. 

Stars #3 and #4 show a Br7 line and a weak Hel 
line at 2.11 /im in absorption. A Hel line at 2.05 /im 
is in emission in spectrum of #3, while in absorption in 
the spectrum of # 4. From a co mparison of the spectra 
with t hose given in lBibbv et al.l (2008) and lHanson et al.l 
(1996), star #3 appears to be a B1-B3 supergiant. Star 
#4 is probably earlier than #3 (09-BO) because it is 
fainter in the Xs-band and the Br7 absorption is weaker. 

The spectra of stars #2, #7, #9, and #11 do not 
show any lines, however, the absence of CO band-head 
suggests a spectral type earlier than G type. 

3.2. Color-magnitude diagrams 
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Fig. 5. — Measurements of EW(CO) versus spectral type for 
the late-type stars identified in GLIMPSE9 (asterisks). For com- 
parison, known giant stars (triangles) and RSG stars (diamonds) 
are also shown dKleinmann fc Hall 1986). The lines indicate a lin- 
ear fit to the EW(CO) versus spectral type of giants (dotted) and 
RSG stars (continuous). There is a small difference between the 
effective temperature of a K5 and an MO star (~150K). Further 
subclassification using the EW(CO) are not possible. 

Color-magnitude diagrams (CMD) of HST/NICMOS 
point sources are presented in Fig. [51 The CMD of the 
cluster field presents a clear sequence of stars with H- 
Ks=^l mag and Ks< 15 mag, while only a few stars 
populate the same region of the CMD for point sources 
extracted from the control field. The control and cluster 
fields were observed in the same way, and therefore have 
identical area. 

Magnitudes in the 2MASS system are preferred for 
the CMD in order to have a dire c t comparison wit h 



the CMDs in 

IDavies et all (|2007D . 



Messineo et"all (|2009t ) . iFiger et all (|2006f l. 



To isolate the cluster sequence, we performed a sta- 
tistical decontamination using stellar counts per 0.5 mag 
bin of (H—Ks) color and 1.0 mag bin of Kg magnitude in 
both cluster and field regions; we randomly subtracted 
from each cluster bin a number of stars equal to that 
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of the corresponding field bin. The resulting diagram is 
shown in the right panel of Fig. [6l 

There is a kink in the "clean" CMD at iCs oc 15 - 16 
mag. This could be due to a poor field subtraction, or 
to a pre-main sequence. If we presume the faint stars 
to be a pre-main sequence, then their ages would range 
between 0.5 and 3 Myr (Fig. O. NGC7419, which con- 
tains 5 RSGs and ha s a age of about 10 Myr, a lso shows 
a younger sequence (jSubramaniam et "aLll2006D . Further 
observations are needed to study the nature of these faint 
stars. 

From the CMD, we estimate an interstellar extinc- 
tion of ^Ka= 1-6 ± 0.2 mag by measuring the median 
H-Ks color of cluster with Ks < 15 mag and 0.6 < 
H—Ks< 1.4 mag, and by adopting the extinction law by 
iMessineo et al.l ()2005[ ). This measurement is indepen- 
dent of age because in the H-Ks versus Ks diagram, the 
isochrones are almost vertical lines. 

We spectroscopically detected several massive stars: 
three RSGs, with spectral type from MO to M2, and two 
blue supergiant stars (BSGs) (09-B2). Since RSG stars 
span a broad range of magnitudes, they cannot be used as 
a distance indicator. Therefore, to determine the cluster 
distance and other dependent parameters (e.g. age and 
mass), one must rely on stars other than RSGs. For the 
BSGs #3 and #4 we measured an interstellar extinction 
Ak,= 1.3 mag and 1.5 mag, respectiv ely (using intrinsic 
magnitudes and exti nction law from iBibbv et al.l 120081 : 
IMessineo et al.l[2005l ). These extinction values are consis- 
tent with that of the cluster sequence, and suggest mem- 
bership. For star ^3 we obtained a spectrophotometric 
distance of 4.2 ±0.4 kpc, and for star #4 of 4.7 ±0.4 kpc. 
These values are consistent with a distance of 4.2 ± 0.3 
kp c, which was i nferre d for the GMC and the SNR W41 
bv iLeahv et al.l ()2008[ ). In the following we will adopt a 
distance of 4.2 kpc. 

The three brightest stars (6.3 < Ks< 7.2 mag) 
have values of EW(CO) typical of RSGs. By compar- 
ing their observe d colors with intrinsic colors of RSGs 
jKoornneef '1983'), and using the extinction law given in 
[Mcssinco ct al. (2005), we estimated the following val- 
ues of interstellar extinction: Ak^— 0.6 ± 0.1 mag for 
star #1 (i^s=7.17 mag), ^k^^ 1-5 ± 0.1 mag for star 
#8 (i^s=6.3mag), and ^k,= 1-7 ± 0.05 mag for star #5 
(iis=7.05 mag). Stars #5 and #8 have extinction val- 
ues consistent with those of the cluster sequence, and are 
likely members, while the bluer color of star #1 suggests 
a foreground star. 

Using the relation between effective t emperature and 
bolom etric correction for RSGs given by iLevesque et al.l 
()2005f l. and a distance of 4.2 kpc, we derived a bolo- 
metric luminosity Mbol=— 5.7 ±0.1 mag for star #8 
and Mbol=— 5.1 ±0.1 mag for star #5. From non- 
rotating evolutionar y tracks with Solar abundance by 
iMevnet et all (|1994[ ). we inferred masses of 10.5 ± 1.5 
Mq, which corresponds to an age of 22.5 ±4.5 Myr. Sim- 
ilar range i s obtained when us i ng the newer non-rotating 
models bv lMevnet fc Maedeil (|2003f ). By increasing the 
distance by 40% (6 kpc) the minimum age would de- 
creases by ~ 30% (15 Myr). 

Star #8 shows water absorption at the blue edge of 
the K-band (1.9-2.1 ^m). Its EW(CO) argues for a late 
giant type (M7III), or an M2 RSG. The water absorp- 
tion indicates the presence of a circumstellar envelope. 
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Fig. 6.— Left panel: HST/NICMOS H-Ks versus Ks color 
magnitude diagram of point sources extracted from the cluster im- 
ages. Middle panel: Similarly, a color magnitude diagram for 
the control field (of equal area) is shown. Right panel: A color- 
magnitude diagram of the cluster field after statistical decontam- 
ination. The ver tical line indicates an iso chrone of 15.0 Myr and 
solar metallicity dLeieune fc Schaerej|200ll) . which was shifted to a 
redd ening of Aks=1.5 mag and a distance of 4.2 kpc (Lcah ~ st al.l 
12008). Dotted lines show isochrones corresponding to a population 
of 0.5 Myr, and 3 Myr with solar metallicity (Siess ct al. 2000). 
The diamond symbols indicate the location of RSGs. Squares in- 
dicate BSGs, and a plus symbol a candidate early type star. A 
triangle shows a giant star, which is likely unrelated to the cluster. 
Stars are numbered as in Table [T] 

Star #8 has a i^s-lS] = 1.5 mag, redder than the 0.42 
and 1.14 mag of star #1 and star #5 (magnitudes at 8 
^m are from the SPITZER/GLIMPSE surve y). Water is 
typically seen in AGB stars and RSG stars (jTsuiil 120001 : 
iBlum et "all l2003f ) . However, the location in the direc- 
tion of a young cluster, and the rarity of bright infrared 
stars, suggests a RSG cluster member. We estimated a 
stellar density of 0.17 ±0.1 bright stars {Ks< 7.5 mag) 
per square arcminute, in the longitude range from 20° to 
30°, and within 0.5° from the Galactic plane; there are 
three such bright stars at the location of GLIMPSE9, 
and one is likely to be a chance alignment by virtue of 
the bluer color. Radial velocities are needed to firmly 
solve the puzzle. 

3.3. Luminosity function 

The cluster and field LFs of stars detected with the 
F160W filter are shown in the top panel of Fig. [71 while 
that of stars detected with the F222M filter are in the 
middle panel of the same figure. Since the transforma- 
tions to the 2MASS system depend on a color term, the 
VEGA photometric system must be used to indepen- 
dently analyze the individual bands. The cluster LF is 
estimated by subtracting the field stars from the observed 
LF, and dividing it by the completeness factor (see the 
Appendix). The cluster LF shows an evident excess of 
stars with [F222M] brighter than « 15 mag. 

MIPSGAL data reveals that the control field has in- 
creased 24/.tm background emission compared to the clus- 
ter field (see Fig.[T]), indicating that it may have a slightly 
higher extinction. However, from the CMDs and LFs it 
appears that this does not affect the analysis for stars 
brighter than [F222M]< 15 mag {Ks<-^ 15 mag). 

We also construct a LF for an extinction free magni- 
tude defined as 
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m-e = Ks - 1.5 X {H - Ks- 1.0), 

where the 1.5 constant is the ratio between interstel- 
lar extinction in j^Ts- band and the reddening in H—Ks 
essmeo etalllMHl) . and 1.0 mag is the average H—Kg 
of the cluster sequence. By using me, each point source 
is moved on the CMD along the reddening vector to an 
observed H—Ks=l mag. Because in H—Kg the intrin- 
sic color of stars i s almost indepen dent of spectral type 
(within 0.3 mag) (iKoornnee l\T98f). the comparison be- 
tween field and cluster luminosity functions becomes ex- 
tinction free. If the control field would have a system- 
atically higher reddening than the cluster field, stars de- 
tected in the control field would have higher H~Ks and 
lower Ks than the true foreground population of the clus- 
ter. This could cause an artificial excess of bright stars 
in the cluster LF of the Kg band, which in turn would 
bias the mass function. The LF with rrie shows a similar 
excess of bright stars (lower panel of Fig. [7|). Thereby, 
the stellar over-density at the location of the GLIMPSE9 
candidate cluster is confirmed. 

A drop appears in the LF around rrie ==12 mag. This 
is due to an excess of stars with {H—Ks)< 1.8 mag in 
the control field, and disappears when applying a color 
selection before building the mass function. 

3.4. Mass function 

Apparent magnitudes of cluster members can be trans- 
formed into initial masses by assuming an interstellar 
extinction, a distance, and using an isochrone. We con- 
sidered Ak^ = 1.6 ± 0.2 mag, as estimat ed from the 
CMD (see above), and a distance of 4.2 kpc (|Leahv et al.l 
120081) . In addition, we used non-rotational models from 
the Geneva group, with increased mass-loss, Solar abun - 
dance, ages of 15, 20, and 27 Myr (jMevnet et al.lll994f ). 
and with the color t r ansfo rmation for Johnson filters by 
iLeieune fc Schaereil (|2001h . The relations between the 
actual mass and the apparent if-band magnitude for a 15, 
20 and 27 Myr p opulations are shown in Fig. [H We used 
the equations bv lKim et al.l (2005) to transform the the- 
oretical isochrones for the F160W, F222M, and Kg short 
filters. An average difference of 0. 07 ± 0.05 mag is f ound 
between the Ks magnitudes from lKim et al.l (|2005l ) and 
those with the empirical transformation given in Sect. 
2.2. This is a measure of the transformation uncertainty. 
Mass functions for the GLIMPSE9 cluster for different 
ages and bands are shown in Figs. [51 and [TUl 

We built mass functions for stars detected in the 
F160W mosaic, as well as for stars in the F222M mosaic. 
We selected a mass range log(M/M0) from 0.1 to 0.8- 
0.9 to ensure that stars were above a completeness limit 
of 80% (see the Appendix), and to include only main- 
sequence stars (see Fig. |8]). Evolved stars fall in a single 
bin because of the degenerate mass-luminosity relation 
(Fig. [HI). For an age of 15 Myr, a linear fit to the mass 
function in the mass range from log(M/M0) = 0.1 to 
0.85 with a bin size of log(M/Mo)=0.05 yielded a slope 
of —0.77 ±0.20 when using the F160W magnitudes, and 
of -0.70 ± 0.25 when using the F222M magnitudes. For 
an age of 20 Myr, the average slope from the two bands is 
-0.63 ± 0.05, while for an age of 27 Myr is -0.64 ± 0.09. 
Bin size variations from log(M/M0)= 0.035 to 0.07 yield 




10 15 
F160W[mag] 




F222M[mQg] 




[mag] 



Fig. 7. — Luminosity functions. Top panel: Observed LFs 
of stars detected in the F160W mosaics. The thick continuous 
line shows the GLIMPSE9 cluster LF after field subtraction, while 
the thick dashed line shows the same LF after correcting for incom- 
pleteness. Similarly, the thin lines show the LFs of the control field 
without (continuum) and with (dashed) correction for incomplete- 
ness. Middle panel: Observed LFs of stars detected in the F222M 
mosaics. Lower panel: LFs with the extinction free magnitude 
mc of stars detected in both bands. Since me depends on color, we 
do not apply any extinction correction. Notice that Glimpse 9 has 
more luminous (massive) members. 

slopes with a standard deviation within 0.1, i.e. the un- 
certainty due to bin size variation is smaller than the 
uncertainty of the fit. 

Besides this uncertainty due to age, the error in 
subtracting the foreground and background population 
needs also to be considered. If the control field is not 
representative of the background and foreground pop- 
ulation seen towards the cluster, a systematic error is 
introduced when subtracting the field stars. The incor- 
rect field subtraction would propagate into the LF and 
mass distribution. As seen in Fig. [TJ the control field 
happens to be located in a dustier region. From the 
CMDs, the field appears to have an excess of red stars 
(2 < H—Ks< 2.4 mag), which implies an extra extinc- 
tion of >lKa=^ 0.5 mag. When artificially reducing the 
attenuation of the control field by j4k= —0.5 mag, the 
slope of the initial mass function decreases by ~ 0.33 
(-1.1 ± 0.1, -1.0 ± 0.16, -1.0 ± 0.1 for 15, 20, and 27 
Myr). We also calculated a mass function with the ex- 
tinction free magnitudes me, which takes into account 
differential extinction, and measured a slope of — 1.0±0.3, 
-0.96 ± 0.24, and -0.92 ± 0.27 for 15, 20, and 27 Myr. 

The origin and universality of the initial mass func- 
tion remain under very active investigation and discus- 
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15 Myr 

20 Myr 
27 Myr 



F222M[mag] 

Fig. 8. — Actual mass versus apparent [F222M], for populations 
of 1 5, 20, and 27 Myr, with Solar abundance and increased mass- 
loss l|Leieune &: Sc haercr 2001). A distance of 4.2 kpc, and ylK=1.6 
mag are used. 

sion. Typically, it is as sumed as a po wer law with an 
exponent of F = —1.35 ()Salpeteilll955f l. Studies of star- 
burst clusters, such as the Arches cluster, r eport flatter 
functions (e.g. -0.7 ±0.1 iFiger et al.lll999D. Recent re- 
determinations of the Arches MF slope bv'Espinoz a et all 
([l009) indicate a slope of —1.1 ± 0.2, which is closer to 
a Salpeter. 

The GLIMPSE9 cluster shows a mass function slightly 
flatter than the Salpeter's one. We tried several ways 
to measure the mass function, including harshly de- 
reddening the control field, and yet we consistently 
come up with slopes that are shallower than Salpeter. 
Mass segregation, which bring massive stars to sink 
into the cluster center and low mass stars to disperse 
into the field, s eems a plausible explanation for this 
()Kim et al.ll2006f ). However, considered the uncertainties 
in age and background subtraction, and the fact that the 
GLIMPSE9 cluster is several million years old, a Salpeter 
initial mass function cannot be excluded. 

3.5. Cluster mass 

By integrating the masses of the candidate member 
stars (using our background-subtracted mass function) 
down to a mass of 1.0 Mq, we measured a cluster mass 
of 1500 ± 300 Mq, where the error takes into account 
the age uncertainty. Since systematic errors could be 
due to errors in the field subtraction, we also estimate 
the cluster mass by integrating only to Ks= 15 mag, 
i.e., by using the upper part of the diagram where field 
contamination is negligible, and extrapolating to 1.0 Mq 
with a power law. When using the slope we estimated 
from the mass distribution, a mass of 1600 ± 400 Mq is 
obtained, while with the Salpeter mass distribution the 
mass value rises to 2100 ± 200 Mq. 



4. CLUSTER SURROUNDINGS 



4.1. 



A giant molecular cloud with SNRs 

The GLIMPSE9 cluster is located at (l,b)=(22.76 ° , 
-0.40 °), in the direction of a GMC (jPame et al.ll2001f ). 
The CO emission from the molecular cloud peaks near 
1=23.3 °, b=-0.3 °; it extends over two degrees in longi- 
tude with a Hnc -of-sig ht velocity from 70 to 85 kms^^ 
(jPame et al.l 120011 ) . lAlbert et al.l (|2006f ) derived an up- 
per limit for the total H2 mass of 2.1 x 10^ M©, assuming 
a (near) kinematic distance of 4.9 kpc. The linear size 




-0.2 0.0 0.2 0.4 0.5 o.a 1.0 1.2 
io5[M;i„H„i;] 





-0,2 0,0 0,2 0,4 0,6 0,8 1-0 1,2 
l05[M(1nllioi;] 




Fig. 9. — Background-subtracted mass functions of the 
GLIMPSE9 cluster from [F222M] or [F160W] magnitudes, for ages 
of 15, 20 and 27 Myr, are shown with a continuous line. The same 
function after correction for incompleteness is shown with a dashed 
line. Error bars represent the statistical Poisson uncertainty. Lin- 
ear fits to bins with log{M/M0) from 0.1 to 0.85-0.90 (i.e., F222M 
between 17.5 and 6.0 mag) are also shown. 

of about 100 pc is in the range of other known giant 
molecular clouds. Several SNRs coincide in projection 
with this molecular clo ud (see Fig. [TT|) . Two SNRs are 
listed in the catalog of lGreenI (fl991l) . G022.7-00.2 and 
G023.8-00.3 (W41). Three other candidate SNRs were 
identified by Hclfandiiin (2006), G22.75 83 - 0.4917, 
G22.9 917 - 0.3583, and G23.5667 - 0.0333. iLeahv et all 
(|2008f) concluded that the SNR W41 is associated with 
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TABLE 2 

Positions of the supernovae from literature. 



ID 


NAME 


RA(J2000) 


DEC(J2000) 


velocity (km s 


diameter (') 


References 


1 


SNR/W41 


18 34 46.42 


-08 44 00 


77±5 


30.0 


Green (1991): Leahy et al. (2008) 


2 


SNR22. 7-0.2 


18 33 17.86 


-09 10 35 




30.0 


Green (1991) 




G022. 8-0.3 


18 33 45.50 


-09 09 17 


82.5 


10.9 


Kuchar & Clark (1997) 


3 


G'22.7583 - 0.4917 


18 34 28.30 


-09 16 00 


74.8 


4.7 


Kuchar fc Clar^ (1997 ) 






18 34 28.00 


-09 16 00 


76.0 




Bronfman et al. ^996) 






18 34 26.70 


-09 15 50 




5.0 


#33 in Helfand et al. (2006) 


4 


G22.9917 - 0.3583 


18 34 26.59 


-09 00 09 




4.5 


#34 in Helfand et al. (2006) 




G22.9 - 0.3 


18 34 12.60 


-09 01 20 


70.9 


4.8 


Kuchar & Clark ( 1997) 


5 


G23.5667 - 0.0333 


18 34 17.09 


-08 20 21 




9.0 


#35 in Helfand et al. (2006) 




G23.5 - 0.0 


18 34 19.60 


-08 22 17 


91.3 


6.1 


Kuchar & Clark (1997) 



Note. — Positions and dimensions were measured on the MAGPIS image. 



-0,3 0,0 0,2 0,4 0,6 0,6 1,0 1,2 
l.g[U(lnlliol)] 





-0,2 0,0 0,2 OA 0,6 0,3 



Fig. 10. — Background-subtracted mass functions of the 
GLIMPSE9 cluster from the extinction corrected nie magnitudes. 
Symbols as in Fig. [9] 



the GMC, and h a s a ra dial velocity vlsr — 77.0 kms ^. 
iKuchar fc ClarkI (|1997f) report vlsr = 82.5 kms^i for 
G022.8-00.3, and ^^..^ b = 91.3 kms^i for G23. 5-0.0. 
[Bronfman et al.l (|1996[ ) measured a vlsr — 76.0 kms~^ 
towards G22. 7583-0. 4917. This suggests that at least 
three candidate SNRs are associated with the same 
GMC. In Table [2] we report positions of the candidate 
SNRs and associated line-of-sight velocity, when avail- 
able. CO observations confirm a molecular cloud with 
a peak at Vlsr = 76 ± 5 kms~^ and a fuU-width-half- 
maximum of 22 kms-i at the location of the GLIMPSE9 
cluster (jPame et al.ll2001[ ). 

The classification o f G022.7— 00.2 is reported as uncer- 
tain bv lGreenI (|1991l ). and the same rad io source is listed 
as an HII region in othe r works (e.g. iKuchar fc ClarkI 
ll997HPaladini et al.ll2003[ ). We use d archival radio d ata 
at 20 and 90 cm from the MAGPIS (White et al. 2005) to 
measure the radio spectral index of the candidate SNRs 
(see Fig. [TT] and Table [3]). Given the negative spectral 
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Fig. 11. — 90 cm images from the M ulti- Array Galactic Plane 
Imaging Survey (MAGPIS) llWhite et a l. 20051. Grey circles show 
the locations of the SNRs listed in Table [2] White squares show 
the position of regions listed in Table \3\ for which we calculated 
a spectral index. Longitude is along the x-axis, while latitude is 
along the y-axis. 

indexes, the radio emission appears dominated by syn- 
chrotron emission in all cases. Since G22.9917 — 0.3583 
and part of the radio shell of G022.7— 00.2 overlaps with 
GLIMPSE 8 /im emission, they are probably a compo- 
sition of SNRs and HII regions. Assuming pure circular 
motion, a peak velocity of 78 kms~^, a distance of 7.6 
kpc for the Galactic Center, and a Solar V=214 kms~^, 
a systematic deviation of 5 kms~\ as we ll as a random 
deviation of 5 km s~^. [Leahy et all (|2008l ) derived a dis- 
tance d=3.9-4.5 kpc for SNR/W41. Using an average 
distance of 4.2 kpc, the angular sizes of G022.7— 00.2 
and G023.3-00.3 (W41) ( -27') yield a linear size of 
^^33 pc, while the angular sizes of G22.7583 — 0.4917, 
G22.9917-0.3583 and G23.5667-0.0333 (about 5 ') yield 
a linear size of ~6pc. These physical sizes are well withi n 
the rang e of other Galactic SNR s (jStupar et al.ll2007f) . 
Recently, iBrunthaler et al.l (|2009| ) measured trigonomet- 
ric parallaxes of 4.6 and 5.9 kpc with two methanol maser 
sources found towards G23.01-0.41 and G23.44-0.18, 
suggesting different complexes arranged along the line of 
sight. 

4.2. Stellar candidate clusters. 
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TABLE 3 

Spectral index from 90 and 20 cm MAGPIS data. 



Region 


RA(J2000) 


DEC{J2000) 


diameter{") 


Spectral index 


1 


18:32:37.114 


-09:13:32.15 


77 


-0.82 


2 


18:33:09.794 


-09:12:45.87 


384 


-1.23 


3 


18:33:20.950 


-09:20:17.81 


154 


-0.98 


4 


18:33:56.335 


-09:09:13.63 


154 


-0.78 


5 


18:34:15.082 


-08:20:41.93 


307 


-0.35 


6 


18:34:20.244 


-08:55:07.61 


154 


-1.44 


7 


18:34:26.988 


-09:00:23.31 


154 


-0.36 


8 


18:34:29.686 


-09:16:03.12 


230 


-0.21 


9 


18:34:38.220 


-08:33:10.46 


154 


-0.72 


10 


18:34:47.170 


-08:43:24.00 


154 


-1.18 



■ m 



-20 



The presence of three SNRs with similar velocities 
{vlsr =75-85 kms~^) suggests that massive star forma- 
tion resulting in the production of several massive O stars 
has been active in multiple sites of this GMC. The stel- 
lar cluster GLIMPSE9 represents one episode of this star 
formation. Two other candidate clusters are reported in 
literature in the direction o f the same mole cular cloud. 
The [SDS'2003]117 cluster ()Bica et al.ll2"00l appears to 
be located onto the SNR shell G22. 7583-0.4917, and the 
GLI MPSEIO cluster ( the candidate number 10 in the list 
by Mercer et al.ll2005[ ) onto the SNR/W41 (see Fig.[T2|. 

We searched for other stellar over-densities in the di- 
rection of this GMC. Detection of stellar over-densities 
are difficult due to the patchiness of the interstellar ex- 
tinction. With the exception of GLIMPSE9 cluster, no 
clear over-densities were detected in the 2MASS im- 
ages. To overcome interstellar extinction, we built a 
density map of point sources detected at 3.6 /im, which 
resulted in many spurious clumps. When looking at 
the SPITZER/GLIMPSE images, however, several neb- 
ular emissions are seen in all four IRAC channels of the 
SPITZER/GLIMPSE survey, indicating the presence of 
HII regions. An increased number of bright 3.6/L(m stars 
also appears associated with some of these regions. We 
therefore visually selected regions of nebular emission 
in all four IRAC channels or apparent over-densities of 
bright stars at 3.6/im. A list of the selected regions is 
given in Table [H 

We show 2MASS CMDs of these selected regions, to- 
gether with a comparison CMD of the GLIMPSE9 clus- 
ter, in Fig. [T21 

Several stellar branches are seen in each CMD, and 
each branch appears broadened by differential redden- 
ing. A bluer sequence with H-Kg—^ 0.4 mag is visi- 
ble in all CMDs. A redder sequence with H-Ks—^ 1 
mag, i.e., similar to that of the GLIMPSE9 cluster, ap- 
pears only at certain locations. Assuming a distance of 
4.2 kpc for the GMC, and considering an average of 1.8 
m ag of visual extinction per kpc, and the extinction law 
bv lMessineo et all ()2005f ). a stellar population associated 
with the GMC must have a minimum interstellar extinc- 
tion of ^K,=0.7 mag. The bluer branch is probably due 
to a young stellar population associated with a closer 
spiral arm. Based on the similarity of interstellar extinc- 
tion, we suggest that branches seen at H-Ks—^ 1 mag 
are due to a stellar population associated with the GMC, 
similar to the GLIMPSE9 cluster. 

Region 1, in the direction of SNR G23.5 — 0.0, shows 
nebular emission, but its 2MASS CMD does not show a 
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Fig. 12. — The greyscale shows a 3.6 /^m GLIMPSE image. 
Circles show the locations of the four SNRs. Squares and labels 
show the regions in Table l4l Longitude is along the x-axis, while 
latitude is along the y-axis. 

stellar sequence with similar color as that of GLIMPSE9. 

Region 2, in the direction of the SNR/W41, shows a 
peak of nebular emission in all four IRAC channels, in- 
dicating an HII region. The 2MASS diagrams also show 
a branch of bright stars at H-Ks=^ 0.8 mag. 

Region 3 shows a concentration of bright stars; how- 
ever, no clear sequences are seen in the CMD. 

Region 4 is also in the direction of SNR/W41. This 
region includes th e GLIMPSEIO cluster, for which 
iMercer et al] (l2005h gives a radius of 0.8'; GLIMPSEIO 
coincides with a nebular peak emission, but the actual 
extension of the emission region has a radius of about 5'. 
Several bright stars are detected in this region, and the 
2MASS CMD shows a sequence at H-Ks=-^ 1.0 mag. 

Region 5 coincides with the area covered by the SNR 
G22. 9917— 0.3583, as seen in the 90 cm image. A se- 
quence at H-Ks=^ 1.0 mag is detected. 

Region 6 is a peak of a nebular emission that extends 
and connects to region 5. The CMD lacks stars in the 
redder sequence. 

Region 7 coincides with the SNR G22. 7583-0.4917. 
This region inclu des the [i3_DS'2003]117 cluster area by 
iBica et al.l (|2003l ). 

Two other regions, region 8 and region 9, were ran- 
domly selected, as a comparison fields. Region 8 does 
not show associated nebular emission, and is in the di- 
rection of the SNR/W41. The CMD shows the blue se- 
quence, but not a clear red sequence. Region 9, which is 
at the outer edge of the complex, does not show nebular 
emission, and its CMD lacks a sequence at H-Ks=^ 1.0 
mag. 

5. SUMMARY 

Investigation with HST/NIC MOS data co nfirms that 
object number 9 in the list bv iMercer et al.l lj2005.) is a 
stellar cluster with a well defined sequence in the H—Kg 
versus Kg diagram. Low-resolution spectroscopic obser- 
vations in K-band yield the spectral types of the brightest 
candidate members, and confirm the presence of massive 
stars. Three RSGs are detected, two of which are candi- 
date cluster members, and two BSGs. A spectrophoto- 
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TABLE 4 
Cluster positions. 



ID 



RA(J2000) DEC(J2000) Radius(') 



References 



GLIMPSE9 


18 


34 


09.59 


-09 


13 


53 


0.3^ 


Mercer et al. (20051 


[BD52003]117 


18 


34 


27.00 


-09 


15 


42 


0.6*^ 


Bica et al. ("2003") 


GLIMPSEIO 


18 


34 


47.00 


-08 


47 


20 


0.8*^ 


Mercer et al. f2005~) 


Regionl 


18 


34 


15.08 


-08 


20 


42 


1.2'= 




Region2 


18 


34 


41.09 


-08 


34 


22 


4.0^= 




Regions 


18 


35 


32.22 


-08 


41 


56 


1.2^= 




Region4 


18 


34 


31.59 


-08 


46 


47 


5.0^= 




Regions 


18 


34 


20.00 


-08 


59 


48 


5.0^= 




Region6 


18 


33 


36.03 


-09 


10 


01 


2.7^= 




Region? 


18 


34 


27.69 


-09 


15 


52 


3.3^= 




Regions 


18 


35 


14.46 


-08 


50 


34 


5.0^= 




RegionQ 


18 


33 


35.88 


-09 


19 


08 


5.0^= 





Note. — Coordinates are followed by radius and references. Region #8 and 
#9 are randomly selected. 

" The radius is measured in the ifg-band image as the half light radius. 

^ The radius is given in the referenced work. 

" The radius is measured in the 3.6 /^m-band image. 
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Fig. 13.— 2MASS J-ft's versus Ks, and U-Ks versus Ks 
diagrams of selected regions in direction of the CMC (see Table 
|4]l. For clarity, an area of 3' is used for all CMDs. Two vertical 
lines at (H—Kg= 0.6 and 1.4 mag indicate the color region of the 
GLIMPSE9 sequence. A similar sequence is visible in reg #2, #4, 
#5 and #7. 



Fig. 13.— [131 continued. 

metric distance of 4.2±0.4 kpc is derived, an age between 
15 and 27 Myr, and a mass of at least 1600 ± 400 Mq. 
The cluster mass function appears slightly flatter than 
the Salpeter's one, which could be the effect of mass seg- 
regation. 

The cluster is located in the direction of a molecu- 
lar complex, which hosts several SNRs and HII regions. 
The cluster distan ce agrees well with that inferred for 
the complex by Lea hy et al.l ()2008[ ). A stellar population 
possibly associated with the giant molecular cloud is seen 
in four other regions: two regions are associated with the 
SNR/W41, one region with the SNR G22.9917-0.3583, 
and another region with SNR G22. 7583-0.4917. The de- 
tection of massive stars in GLIMPSE9, and the concomi- 
tant presence of several SNRs render this GMC particu- 
lar interesting. It is a good laboratory to investigate var- 
ious issues about massive stars and multi-seeded star for- 
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mation. By identifying stellar clusters of the same com- 
plex, one can study their properties and variations across 
space and time. Some of the most massive stars of the 
Milky Way may be hiding among the cluster members. 
The detection of short lived massive stars (e.g. Wolf- 
Rayets, Red Supergiants, Luminous Blue Variables) is 
important to understand their formation, evolution and 
fate. By association with the SNRs, massive stars also 
yield the initial masses of the supernova progenitors. A 
more detailed study of this complex will be presented in 
a future paper. 
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APPENDIX 
ARTIFICIAL-STAR EXPERIMENTS 

For a proper cluster analysis (e.g. for a statistical field decontamination, luminosity function, and mass function 
studies), one needs to characterize all undesired biases of photometric measurements in dense stellar fields. Crowding 
limits the detection of point sources, and such incompleteness varies from field to field because it depends on the stellar 
density. We, therefore, run simulations with artificial stars. 

Artificial stars were created with the observed PSF. Then, the artificial stars were added into each mosaic at random 
positions, but imposing only one star over a 50 x 50 pixels area in order not to alter the level of crowding. The 
luminosity function of the artificial stars was assumed identical to the observed luminosity function of each frame. The 
photometry of the artificial stars was recovered following the exact steps as those performed for the actual catalog. 
The procedure was iterated 300 times, giving a total of 27,000 artificial stars per mosaic. By comparing the fluxes of 
the artificial stars' input with their fluxes after re-extraction, we obtained an estimate of the photometric uncertainty 
and completeness limit. 

In Figure [Ml the differences between the input (maginput) and the extracted magnitudes (magoutput) of the artificial 
stars are shown as a function of maginput] differences remain, on average, null till a certain magnitude, then increase 
with increasing maginput', this is due to blending between artificial stars and real stars, and it causes a bin-to-bin 
migration in the luminosity function (LF). A star is considered lost if no star is located by the DAOPHOT point- 
source finding algorithm within 1.5 pixels (0.1") of the inserted location. The fraction of recovered stars per bin of 
magnitude is the completeness factor (C/). A C/ above 80% is reached for [i^222M] < 17.0 mag in the cluster field, 
and [F222M] < 18.0 mag in the control field. For the cluster field, the one sigma deviation is below 0.3 mag for 
[i^222M] <^ 16.0 mag, but 1.0 mag for [F222M] =~ 18.5 mag. In the control field, the one sigma deviation is below 
0.3 mag for [F222M] < 18 mag, and 1.0 mag for [i^222M] =~ 20 mag. 
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